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SUMMARY 


An analysis of a solid propellant ranjet which carries the 
propellant in the contmstion chamber was made, and a procedure 
outlined for calculating the required burning rate of such a ranjet 
which develops constant thrust. The two factors which influenced 


this develovnent were: 


1) The solid propellant ramjet's combustion chamber varies 
in size with time due to burning of tho fuel. 


2) The solid propellent ramjet develops constent thrust. 


These factors wore coordinated with the analysis of the internal 
flow system of the ramjot, ond an expression for the required 
buming rate was derivod. 
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skin friction coefficient 


specific heat at constant pressure 
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I. INTRODUCTION 


The use cf solid fuel in a ramjet introduces a new problem, 

Ag she fuol burns, ths effective combustion chamber area a ا‎ 
This would normally result ኤጩ change in thrust; however, if the 
concition cf constant thrust is imcosed, then this change in 
com™vustion chamber area would wary the amount of heat required. 

ihe chance in heat required may be expressed in terms of the burning 
rate required for a soiid provellent reanjet developing constant 
thrust. 

References 1 and 2, which are representative of the numerous 
analyses in this field, were concemoed with the intemal flow throuch 
a ranjet having = constent ccnbussicn chamber ares, and no restriction 
on the thrust. These analyses were developed for subsonic conditions, 
and a limited heat addition. 

In the present investigation, a theoretical amslysis was mde 
for supersonic conditions of a solid fuel ramjet, which carries the 
vrorelians% in the combustion chanber, developing constant thrust. 
This latter requirenent fixes the exit conditions and permits 
calculation of the intemal flow from the rear of the ramjet to the 
combusticn chamber inlet. in this manner, the pressure icss 
asross tha combustion chamber and nozzle was computed. The pressure 
risg across the diffuser may be determined from the mown inlet 
conditions. By equating the rressure rise to the pressure loss, 
uniane determination of the heat required for constant thrust can be 


made. na amount of heat developed 4s a function of the fuel and 





سا۔ےہ 


its distribution, end from this, th: rewired burning rate may 
ba determined. 

osulte of this ۳919 وحه‎ presented in Dart II of the 
investivotion in tha form of worlsine sgu2ztions, some of which have 
Loon solved end presonted on chorta, These equations and charts 
haya vocon consolidated, and tha romíred durnine rate or specific 
fusl consumption is presented in Figure 10. Pert III illustrates 


the method of solving a numerical problen. 
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Figure 1 Schematic of a solid propellant ramjet which 
carries its fuel in the combustion chanber. 





IJ. ANALYSIS OF INTERNAL ILOW SYSTEMS 


Flow cheracteristies of the intemal flow system must be 
calculated so that the state of air entering the combustion chamber, 
the flow change ecross the combustion chamber, the net drag of the 
internms1 flor, and the required geometry of the ramjet can be 
determined. In determining the intemal flow from the rear end of 
the ranjet, the following relationshins mst be esteblisheds 
(A) nozzle throat condition to ambient condition; (B) flow changes 
across the nozzle; and (C) flow changes across the combustion 
chamber. These expressions must then be rolated to (D) flow chances 
across the diffuser. 

One dimensional flow, uniform velocity distribution across the 
duct, constant ees properties, end idesl expansion were assumed. 

Initial quantities assumed to be known are as follows: 

(1) The desired thrust at designed altitude. 

(2) The heeting value ond density of the fuel. 


(3) Tesimed Mach number (surersenic) snd atmospheric 
conditions Pos Por To’ 


(4) Skin friction coefficient. 
(5) Aroa of the throat of the nozzle. 
(5) Ratio of length of chamber to throat diameter. 


(7) The number of grains (N ) of fuel. 





A complete investigatio. siiovld include consideration of: 
1) isentropic diffusion, and 2) normal shock wave; and whether 
the nozzle 4s sonic or not at the throat. This analysis cutlines 


in detail the case of isentropic d&ffusion with a sonic nozzle. 


— 209 Relation Down iogsla 244 and Ambient Condi tions 
An exact relation between nogzle exit condition and thet of 
che anbient surroundings may be obtrined from the expression for 
thrust. The thrust of a ramjet is given by the increase in 


nonantun of the internal flow, or 


Yor UVa 


The mass rate of flow at the nozzle throat is 


n= Ps Ve SZ Pe Sy +, since ii, = 1 


tho thrust may then be axorossed as 
D = - 
Pe Py a dA, (v, Ye) 


Hz carl aT, -M, fear) 
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All of the factors on the left hend side of the equation aro 


known initially, therefore a new parameter, E , is introduced 
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کہ‎ 
08 Pe 2 (3 ) 
R R ከ ከ 
end Stil 
[ie TT 
ጌ Tr = ሙሮ 
Using the isentropic relation 
Te _ 7 
Tr It ME 
then ls 
-3(7-9) 
| sr e 
۸۰ = Mo 14 
= 2-۱ 2 
is 
Solving for  ر‎ 
stl 2 
xt! or 
To = Ma - E ጅሮ El My 
Te x Mo 


(LA) 


IA is nob a known parameter, but it can be related to the tomperature 
rise. 


Tho hoat liberated is egual to the total temperature rise, or 





‘Ise - Tat a H 
T Cp lat (2) 





The compression fron free stream to station 2 at the combustion 
chamber inlet is an isentropic one, and the state 2 is assumed to 
be a stegnation state, where 
ہے‎ := ረ 

Similarly the expansion from state 3 to stato 4 is an isentropic ons 
end 

Tye‏ = وو] سے 
Equation (2) becomes‏ 


Tec 


Tor 


=y 


6 
ን. (3) 


Using the isentropic relation 
“fat = La (i+ 5.- Má) 


equation (3) becones 





s-! = 
Be I IIT 
Ta b+ SE |۱۱ + لت‎ (4) 
AH 
Zquating (1A) and (4), and solving for aa » results in the 
p foe 
following expression, 
2: a X (5) 
AH Mo ZE My 
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Tmi3 equ ticn can 
he solvod in tho fom 
of e graph and would be 


ciníler to that shown 
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Figure 2 


for a ronjet with no 


ac 


The ideal value of heat required መ.- ۷ q 
icsses can bo obteined by entering Figure 2 with M, = He As losses 


develop, the oxit mach namber Mı Gecrenses, ond a larger value of 


AH 


= is roquired to maintain constant thrust. 
p tor 


Be Flow Changes Across Tus Nogsle 


In development of the analysis from the rear of the ranjet, the 
flow across the nozzle must be expressed in terms of the combustion 


chembar outlet. The pressure ratio across the nozzle is 


x 
E ۳ 
۴ PR - ) (6) 
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By the continuity eqiotion 








Ma = Me 
Nr DE ٠ Td” 
ዬ a, Ms + 5 (x de Me “ን 
or 
Ei 72 
P aT E مت‎ 
||. = 3 As (TS) 
ALSO, 





T 
Solving for Te ۰ 
2 (5-1) 


= ጀተ | 
13 (Er 4፳%ዝ (8) 
[+ St! Ms Iv De 
Substitution of this relation into equation (6) yields 
Gr | 2 ۳ EN (9) 
R 5. yti My رک‎ 


A 
A relation between wi and M, my be obtained by equating (7) 


ሸተ! 
d. ۳ | 1. 
Ind E 











and (8) 








(10) 





سوت 


In equetien (9), the term a is a known quantity end 
hes physical sicnificance, whoroas Mz 1s an unknown value and is 
rot gmerally used as a parameter. This analysis will therofore 
vo Ceveloped in terns of ከ ; however, because of case of 
conversion, as-seen in equation (10), solutions to all equations 
were derived in terms of Mz. The final result is expressed in the 
known paranoter, 5 . 


äquation (9) may be written as 


X 


Ps ከፍ | (33) 


Im 








The flow through the combustion chamber is considered next. he 
channel or duct has a constant cross section in which friction losses 
occur. As the fuel burns uniformly, the cross section of the duct 
increases, end a change in heat added is required to maintain constant 
tniuste Ii keeps with tho proceónre outlinea, all relations will 
be stated in terms of combustion chamber outlet conditions, 

The pressure drop between stations 2 and 3 (at combustion 
chamber inlet and outlet) is obtained ty applying the momentum relation 


petwoon thease stations: 


۰ De 
ር a n = ہش‎ 
p ے0‎ VAR Pa + P, ሃን ተ A 3 





where Dp is the drag force due to friction within the conbustion 
chamber, and An = An. 

Fron equation (12) it 1s evident that the pressure drop is 
composed of a component associated with the frictionel drag force, 
ent a component dus to the momentum increase in the duct. ^S 
rointed out in Reference (1), the former comronent is frequently 
nezlected but becomes very important when the velocity in the 
duct is large. ‘The drag term is presented in the fom of a drag 


coefficient, Cor , where 
5 6 
Cor 5 dı Az (13) 


iz the incompressible flow is taken as a basis for the unheated low 


speed condition, the drag coefficient becomes 


Ds 
Cori = ( A, Az ,, (34) 


and is determined by the geometry of the ramjet snd the Reynolds number. 
There is en increase in the dynamic pressure when heat is applied, 
resulting in an increase in the incompressible dreg coefficient. 

if the dynamic pressure variea linearly with the length of the 
conbustion chamber, an average value may be assumed, and equation (14) 


becones 
Dei 
ር 41 ይ 651). 


if the incompressible drag is now expressed in terms of the skin 
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friction coefficiont (a known value), Cp; becomes 


„Va ۷ 
ብዋ Lee |(# DS 


EDS ٦ 
4. 2 ۷۵ + (3 V3 Nr Dc 
SE 
or 
ہت‎ = 4 “De Cf 
hence 
گے‎ ጴህ ^ 4V3 
DS A (15) 


Since the heating effect has bean accounted for in the basic drag 


relation, the two drag equations (13) end (15) may be equated 


cn EM A = HS EO + D4, 


or 
Cpg = Da Ge (1+ (s/n, ) (16) 


With this definition of the drag coefficient, equation (12) becomes 





Bfn) 


,3 2 + وا + وم = Ou‏ 


or 





— ے‎ | + rm 522 (1+ o, ) + (| 





-1% 
or, in teras of tho dianoter ratio, since D, v:ries, 


Pr 2 
ها‎ am 


In order to evaluate the pressure ratio, ths density ratio 
must first be obtained. This ratio may be determined by sinulteneous 


solution of tha energy equation, 


2. 2 
Va 


= ۷ en 
888820 


the continuity equation, 


ቴ እ = (3 V3 


the perfect gas relation, 


= 1 
e R 
and the momentum relation, equation (12). The final result is 
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a 
or in terms of 25: * 
۳ ው ር کا‎ | EE Y ی‎ Y 
o aU ay 7) +۲۸۷ OA | + a 


(18) 





Solution of equation (17) provides 2 , for, 





-13=- 


ደየ‏ رھ 
ዴ ٩ P, (19)‏ 





Thus the pressure drop across the combustion chamber and nozzle has 
been calculated, and must equal the pressure rise across the 
diffuser. 

Solution of equation (17) also permits determination of Mz , 
for by the continuity ecuation, 

| 0. a, Ma = fa as Mz 
or 

mo = E = 


Also, 





h ۳ ee 2) = ۳ :.. E 
۳ at + 5: ay Y + کے‎ E 


Conbining the last two equations, and simplifying, 


۳ | 
زا‎ A 
ںہ‎ %4) - ነጻ 


The pressure rise across the diffuser as computed from the rear 
of the ramjet vas plotted vs the conbustion chamber inlet mach 


number as calculated in equation (20), (ses Figure 3). 
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M VS. حت‎ 
= Po 


L 
for constant E, 3156 a 
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Figure 3 


A plot of combustion chamber inlet mach number and heat 


required for constent thrust as developed in equations (18) 


and (20) 49 shown in Meure 4. 


INL 
for constant E, Mo» and TT 





€ plo t 


Figure 4 
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ie Flow Changes Across The Diffuser 


In an anılysis of a ranjet having a designed supersonic mach 
number, there are two types of diffusion to be considered - namely, 
isentropic nå normal shock wave diffusion. In this investigation, 
diffusion was considered to be isentropic, and the diffuser 100% 


efficient. The pressure rise scross the diffuser is 








b E, | EN ME: v- 

3. Ta ig ማመ ۹ 

T) = | = A (21)‏ = سے 
E Mi‏ + | . 


Nou, the pressure rise across the diffuser as computed from the 
front end of the ramjet vas plotted in Figure 5. 
ا ہ۷‎ 


Mi ys: id for constant E, Mo» Ea 





ME 
Pr 
Po 
Figure 5 
s Pp 
Unique determination of Ma. and ے۔‎ for constant 


Ps 
thrust vas obtained by superimposing Figure 5 on Figure 3 and 





--ጊጮ 


notin: ths points of intersection as shown in Figure 6. 


M VS. کے‎ for constant E, Mo, 


Front نا‎ 


Maz 


Rear End 


Figure 6 


A Detemination of Bumine Rate 

The condition of constant thrust has been satisfied by the 
gimultaneous solution of equations (19), (20), and (21). With the 
value of M2 obtained at the intersection of the two curves as 
a parameter, the v>lus of heat required (AE | was 8 
from Figure 4 for a given value of on . After obtaining valuss 
of ne for various en ¢ a chart was next constructed of 


A AH 
era ir (sco Figure 7). 


An n NL 
CT ሃ5- 5ኳቧ for constant E, Mo» ge 
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Tris chart was used later to obtein tho heat required for en 
slonentai volume. In considarinsz an elemental volume, the 
waisht cf fuel consumed and tho time of consunpticn mst be 
calculated, 


To aeeomplish this, the mass rats of flow is written as 


Wir De cere 
mM < 3: ۳۱۸۰ < ۳ ~ a 





(2 &o 
۱7۱ و و‎ 
vıl 
2(r-1) 
202-0 e 6. اچ‎ Mo 
- جم‎ (+ lM 
Then, 
vil 
2(z- 1) 
ሃጥ De 5 t as (22) 
m = pao Ma he d 1+ St my 


The fueleair ratio may be determined from 
AH 


og 


3 
T 


The A A vas obtained from Figure 6, and h is a giver value, 


therefore 
Mme = M 1 


The finel relaticn between the weight of fuel burned, the tins 


of bumins, anù the chanses in diameter may be obtained by 





کے 


considering a smell elemental volume of fuel. Whe volume ( V, ) 


at time T 248 


۷ < De (4 d.) = 


where A De is the change in fuel diameter. At time t+ 4At, 


۷ 
Do Cea በ Det f= ۳ (t) + 1۸۲ (ቲነ (24) 


۷۸۷ De 


Or, in terms of ER ý 





N > (trat በቦ 9+9 NV ኤፌ 
28 رت‎ = a t inr (d 


(N Dc(t) ( ۳۳ ے‎ ) 
= TTT (5 LNT DX (t) 


Hence, solving for Vi, 





ጌ- u en 
۷ = ENT DC (4) Yy de (t) 5 | 








The weight of the fuel consumed, We , mey be obteined fron 


E Vi 





The tins of burning of this elcnental volume is 


We 
ረ! ኮካ ር 





A plot of the woight of fuel consumed vs. diameter ratio and 


vs. time of burning was made in Figures 8 and 9. 


A 





We vs 


In سا‎ we vs. tu for given ሠራ 


Ke Cov quen SE 


a 


Figure 8 Figure 9 


To show the curve in detail, the specific fuel consumption 
minus the ideal specific fuel consumption was plotted. The specific 


fusl consumtion is 


3600 f 3400 ME 


S = TE 5 (26) 


[ 


For the ideal case in which there are no losses, 








_ 2(Y- 1) 
¥ Pp, Ar À ۱ My 
m = 2 MN e vr (27) 
Sp Ter = 





e ÁH 
12 value of Cater sa Obtained from Yiewre 2 entering with 


My œ Mo as a parameter. Then 
3600 AH m 4 
Das አሯ 


where m is as exoressed in equation (27). 
A plot was then constructed of (S- S.) vs. time of bumins 
(596 Figure 10). 
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Figure 10 
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In a namerical problem to determine the buiming rate of a 
solid fuol ranjot (which carries its fuel in the combustion 


chambor) developing constant thrust, tho following deta is given: 


My = 2 h = 14,220.65 BTU/lb. N = 25 

2 = 3 fe = 140.5 2۵۴ Cp = 0.003 
vos 1.4 hy = 1 tt? መ= = 50 

Po = 4.362 pela. و‎ ۸ F = 1758 lb. 


Equation (5) has been solved for these values, and is shown in 
Figure 11. From this chart, the ideal value of heat required for a 


ranjet with no losses is 1.019. As the losses develop, the exit 


AH 


nach number is decreased and a larger value of re is 





required. 

ihe relation between መመ and Ma as expressed in 
equation (10) is shown in Figure 12. 

Solutions of equations (19) and (20) for various values of ma 
aro shown in Figures 134, 138, and 130. It is to be recalled that 
$hese solutions were calculated from the rear of the ranjet. 
Bauation (21), calculated from the front end, is plotted on these 


same firures, and the intersection of these two curves 0 


is 


a má Ma for a constant thrust, solid 


determine 


propellant ranjet. 





-33> 


The heat required to maintain this constant thrust was 
plotted vs. the combustion chamber inlet mach number, end is 
show in Figures 144, 143, and 140. Using the uniquely determined 
vaitts of Ma, as the entering parameter, the heat required was 
obtained as a function of the diameter ratio (E) » and is 
shown in Meure 15. 

Using the process of an elemental volume as doscribed in 
Section ود‎ end solving equation (25), a relation between the weicht 
of the fuel consumed and the diameter ratio was determined. ‘this 
relation is shown in Figure 16. ‘the weight of fuel consumed was 
plotted vs. the burning tîne in Figure 17. The slope of this curve 
is tho required burning rate. Unfortunately, for the values chosen, 
the burning rate is relatively constant. Actually, the required 
burning rate should decrease with an increase in burning time. To 
show the relation between the weight of fuel consumed and the 
burning time in greater deteil, the specific fuel consumption 
was calculated from equation (26) and the value of (3 - (مة‎ vas 
plotted vs. time of burning and ig shown in Figure 18. This curve 
snows a chanse of approximately fifty percent in (S =- So) over tho 


ranse of burning tine. 
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